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Abstract. Complete flux-calibrated spectra covering the spectral range from 6OOOA to 2.5^m have been obtained 
along the HHl jet and analysed in order to explore the potential of a combined optical/near-IR diagnostic applied 
to jets from young stellar objects. The main physical parameters (visual extinction, electron temperature and 
density, ionization fraction and total density) have been derived along the jet using various diagnostic line ratios. 
This multi-line analysis shows, in each spatially unresolved knot, the presence of zones at different excitation 
conditions, as expected from the cooling layers behind a shock front. In particular, a density stratification in the 
jet is evident from ratios of various lines of different critical density. We measure electron densities in the range 
610^-310^ cm"^ with the [Sll] optical doublet lines, 410^-10* cm"^ with the near-IR [Fen] lines, and lO'^-lO'^ 
cm""^ with optical [Fen] and Call lines. The electron temperature also shows variations, with values between 
8000-11000 K derived from optical/near-IR [Fen] lines and 11000-20000 K from a combined diagnostic employing 
optical [O l] and [N n] lines. Thus [Fe n] lines originate in a cooling layer located at larger distances from the shock 
front than that generating the optical lines, where the compression is higher and the temperature is declining. 
The derived parameters were used to measure the mass flux along the jet, adopting different procedures, the 
advantages and limitations of which are discussed. The [Fen] 1.64/jm line luminosity turns out to be more suitable 
to measure Mjst than the optical lines, since it samples a fraction of the total mass flowing through a knot larger 
than the [O l] or [S n] lines. Mjet is high in the initial part of the flow (~2.2 1Q~^ Mq yr~^) but decreases by about 
an order of magnitude further out. Conversely, the mass flux associated with the warm molecular material is low, 
Mh2~10~^ M0yr~^, and does not show appreciable variations along the jet. We suggest that part of the mass 
flux in the external regions is not revealed in optical and IR lines because it is associated with a colder atomic 
component, which may be traced by the far-IR [O l]63pim line. 

Finally, we flnd that the gas-phase abundance of refractory species, such as Fe, C, Ca, and Ni, is lower than the 
solar value, with the lowest values (between 10 and 30% of solar) derived in the inner and densest regions. This 
suggests a significant fraction of dust grains may still be present in the jet beam, imposing constraints on the 
efficiency of grain destruction by multiple low- velocity shock events. 
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1. Introduction 

Highly-coUimated jets of matter from young stellar ob- 
jects (YSOs) are a common phenomena in star forma- 
tion. These flows exhibit rich emission-line spectra over 
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a wide range of wavelengths, from the UV to the radio. 
Many of the important tracers of physical conditions in 
the jet lie in the optical and near IR (NIR). In the op- 
tical, prominent lines such as [Oi]A6300, Ha, [Nn]A6583, 
[S ii]AA6716,6731, have been widely used to derive elec- 
tron density, temperature, and ionization fraction along 
the jet beams (e.g. Bacciotti & Eisloffel 1999, hereafter 
BE99, Hartigan et al. 1994). The NIR spectrum, on the 
other hand, comprises both important atomic lines such as 
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[Fell], Pa/?, [Sii], and [Nl], giving complementary infor- 
mation on the excitation conditions (Nisini et al. 2002), as 
well as intense H2 transitions from various vibrational lev- 
els, tracing the colder molecular part of the post-shocked 
gas (e.g. Giannini et al. 2004, Eisloffel et al. 2000). So far, 
the optical and NIR spectra of the same object have never 
been compared through combined optical/NIR diagnos- 
tics. Moreover, studies in both the optical and the infrared 
have mostly concentrated on limited spectral ranges, e.g. 
around 6500A and 2/xm. Hence the diagnostic capabili- 
ties of lines in other parts of the spectrum have not been 
fully exploited. Only recently has the potential of transi- 
tions from 0.7 to 1.5/xm been highlighted by Giannini et 
al. (2002) and Hartigan et al. (2004). 

The aim of the work described in this paper is to 
fill this gap and build a set of tools for a combined op- 
tical/NIR analysis of spectra of stellar jets. The advan- 
tage of such a multi-wavelength approach is twofold. On 
the one hand it allows us to employ ratios between lines 
from the same species that are well separated in wave- 
length thereby providing more stringent constraints on 
excitation conditions. This is the case, for example, for 
the [Feii]8620A/1.64/im and the [Sii]6716,673lA/1.03/im 
ratios, which are very sensitive to the gas electron temper- 
ature (Nisini et al. 2002, Pesenti et al. 2002). On the other 
hand, such an approach gives us the possibility of probing 
the different components of the jet cooling layers, where 
strong gradients are expected along the jet axis, that are 
not spatially resolved with current instrumentation. This, 
in turn, allows a better determination of parameters, such 
as the mass flux in the jet, which are fundamental to its 
dynamics and other properties. A correct derivation of 
this quantity is particularly important in heavily embed- 
ded sources, where the mass flux is the primary indirect 
manifestation of protostellar accretion activity. 

In addition, a wide wavelength coverage gives us the 
possibiliy of observing lines of less abundant species that 
nevertheless contain important diagnostic information. In 
particular, lines from many refractory species (Ca, Ni, Cr 
in addition to Fe and C) can be used to derive the amount 
of gas depletion in the jet, thus setting constraints to 
the degree of dust grain destruction by shocks. The de- 
gree of grain destruction in shocks, and the relative gas- 
phase abundances among different species, depend on the 
shock properties, but also on the grain structure (Jones 
2000). The latter, however, is not completely understood, 
although the grain composition is of crucial importance in 
many processes, spanning from the physics of molecular 
clouds to the formation of planets. 

For our purpose we have observed the HHl jet by ob- 
taining flux calibrated spectra covering the wavelength 
range from 0.65 to 2.5/um, as described in Section 2. This 
jet, emanating from the radio source VLAl (Strom et al. 
1985) and located in the Orion L1641 molecular (D=460 
pc), is well suited to test combined optical/IR analysis. 
It powers the prototype HH objects HHl/2 and is bright 
both in optical (Eisloffel et al. 1994, Hester et al. 1998, 
Reipurth et al. 2000) NIR H2 and [Feii] lines (Davis et 



al. 1994, Noriega-Crespo et al. 1994, Davis et al. 2000). 
A limited set of optical forbidden, as well as H2, emis- 
sion lines were analysed in the past by various authors 
(Solf & Bohm 1991, Noriega-Crespo et al. 1997, Eisloffel 
et al. 2000, Medves 2003). Here we use our new combined 
optical/near-IR observations to derive the relevant physi- 
cal parameters (Ay, rie, T^, Xe, and nn) and their varia- 
tion along the HHl jet, using different line ratios capable 
of diagnosing regions with different excitation conditions. 
This is described in Section 3, while the diagnostic capa- 
bilities of specific lines are illustrated in some detail in 
Appendix A. In Appendix B we discuss the limitations of 
the widely used [Fell] 1.64/1.25/xm ratio to measure ex- 
tinction along the line of sight. 

In Section 4 we estimate the jet mass flux rate by 
adopting different lines and methods, discussing the limi- 
tations and advantages of the different procedures (Section 
4). Finally, in Section 5, we will set constraints on the effi- 
ciency of shocks to disrupt dust grains by deriving the gas- 
phase abundance of different refractory species. In Section 
6 we summarize our results discussing their implications 
for the understanding of HH jet properties. 



2. Observations and results 

2.1. ESO 3.6m-EFOSC2 and NTT-SofI observations 

Since the aim of this project is to apply a combined op- 
tical/NIR line analysis, we obtained spectra that are as 
homogeneous as possible in the two wavelength ranges. 
Observations with the ESO 3.6-m and NTT telescopes on 
La Silla, Chile, were carried out during adjacent nights 
(7-8 January and 11-12 January 2003 for the optical and 
IR observations, respectively). The instruments, EF0SC2 
for the optical and Sofl for the infrared, had the same slit 
width (1"), corresponding to a similar spectral resolution 
between 500 and 600. The spatial scale of the two cam- 
eras are also comparable, being O'.'314/pixel for EF0SC2 
and O'.'29/pixel for Sofl. The slits were aligned along the 
HHl jet at a P.A. of -36°. The EF0SC2 observations cov- 
ered the spectral rang e from 6015 to 10320A in a single 
grism setting, while the NTT observations were obtained 
with the blue (0.95-1.64y^m) and red (1.53-2.52/im) grisms. 
Integration times were 1800s for EF0SC2 and 2400s for 
each of the two Sofl grisms. Data reduction was performed 
using standard IRAF tasks. In both the optical and IR 
spectra we corrected for the atmospheric spectral response 
by dividing the object spectra by the spectrum of an 08 
star showing telluric lines. The spectra were wavelength 
calibrated with a xenon-argon lamp, giving an uncertainty 
smaller than the spectral resolution element (~ 20 A). Flux 
calibrations were performed using photometric standards 
observed at an air mass difference not larger than 0.3 with 
respect to our target observations. 
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2.2. Extraction of the spectra and intercalibration 

At the end of the reduction procedure, we had three 
individually calibrated spectral images that have to be 
combined to produce single spectra of the same regions. 
The spatial zero point was set to the position of the jet 
source, VLAl. To determine this in our frames, we used 
the bright Cohen-Schwartz (CS) star, adopting for it an 
angular separation of 35'.'9 from VLAl (Rodriguez et al. 
2000). Figurenreports the emission profiles of bright lines 
observed in the three different images, namely [S ii] 673lA, 
[Fell] 1.64/im and H2 2.12 fj,m, that trace different exita- 
tion conditions and/or extinction. From the [Fell] pro- 
files, nine contiguous emission knots have been identified, 
named A to I following the nomenclature of Eisloffel et 
al. (1994) and Bally et al. (2002) 1. Their spatial extent 
is reported in Table 2. We then extracted the spectra of 
the individual knots and computed the line fluxes by fit- 
ting the profiles with single or, where blending is present, 
double Gaussian profiles. Errors to these line fluxes are 
measured from the rms of the adjacent continuum. 

To intercalibrate fluxes for the three spectral segments 
of the same knot, we used the spectrum from the Sofl- 
blue grism as a reference, normalizing the EF0SC2 spec- 
trum to it with the [C ll]AA9824,9850 hnes, and the Sofl- 
red grism spectrum with the [Fell] 1.64^m line. We note 
that the intercalibration factors may differ from unity and 
change from knot to knot, for a variety of causes. These 
include different atmospheric conditions during the obser- 
vations, a small misalignment of the optical and IR slits, 
and slighty different spatial sampling of the instruments. 
The net result were intercalibration factors ranging be- 
tween 0.6 and 1.1 for EF0SC2 and the SofI blue grism, 
and between 0.9 and 1.5 for the blue and red SofI grisms. 
For those knots where the intercalibration lines were not 
observed, average factors, derived for the other knots, were 
used. 

Figure 12 shows the complete, intercalibrated spectrum 
of knot G - the brightest in the optical -and in Table 1 we 
report fluxes and identification for all the detected lines in 
this knot. These, as a whole, indicate low excitation con- 
ditions, as expected for HH objects. In addition, H2 lines 
from the 1-0 and 2-1 vibrational series were also observed. 

2.3. Origin of the observed lines 

The observed spectral range includes transitions from neu- 
tral and singly ionized atomic species which have different 
excitation temperatures and critical densities, thus allow- 
ing us to probe regions with different conditions. In HH 
jets the lines are believed to be excited in the cooling 
region behind a shock front. This zone, which is gener- 
ally about 10^^ cm wide (i.e. O'.'OIS at 460pc) is not spa- 

^ Bally et al. (2002) named the knots in the jet from Aj 
to Ij in order of decreasing distance from VLAl. They also 
found two additional knots, A and B, close to knot Gj but at 
a different position angle, and thus probably belonging to a 
different jet. 



HHl jet 




5 10 15 20 25 



Distance from VLAl (arcsec) 

Fig. 1. Profiles of the [S 11] , [Fe 11] and H2 emission along 
the HHl jet. Vertical dashed lines indicate the range coad- 
ded to obtain individual knot fluxes. The [Fell] profile was 
used to define the knot separation. 

tially resolved in our spectra, but the various lines may 
originate from different parts of the post-shocked regions. 
This is shown in Fig. |3| where we plot the emission pro- 
files (normalized to their peak value) in the cooling re- 
gion of a 70 kms^^ shock. The profiles of Xe, T^, and 
Tig used to compute the intensity profiles were taken from 
the shock model of Hartigan et al. (1994) (their Fig. 1). 
Similar plots were shown by BE99, who noted that the 
optical [Sii], [Nil], and [Oi] transitions originate from 
similar regions at intermediate temperatures and ioniza- 
tion fractions. From our Fig.|21we see that [S 11] transitions 
at about 1.03/im come from a more compact and hot re- 
gion, as they have a higher excitation temperature with 
respect to the optical lines. On the other hand, the emis- 
sion zone of the [Feii]1.64/xm line is broader and covers 
most of the post-shocked region, tracing mainly the gas 
at higher density and lower temperature. The profiles of 
other strong lines, such as [C i]9850A and [Caii]729lA are 
also indicated. In Appendix A, we describe in more detail 
the diagnostic capabilities of the various lines used for our 
analysis. 

In addition to the atomic lines, we also detect H2 emis- 
sion, which is commonly observed to be spatially associ- 
ated with the atomic gas in jets from young stars. The 
H2 lines in the HHl jet all come from the first two vibra- 
tional levels and thus trace temperatures of the order of 
2000 K. The non-detection of lines from higher vibrational 
transitions, together with the relatively low reddening in 
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Fig. 2. Combined EF0SC2 + Sofl spectrum from 0.6 to 2.5/xm of knot G of the HHl jet. The stronger hnes are 
identified. 



at least the external sections of the jet, allows us to in- 
fer that molecular gas at higher temperatures contribute 
little to the overall flux (Giannini et al. 2004). The ori- 
gin of the molecular emission, closely associated with gas 
of much higher excitation, is still not clear. The emission 
could originate in the post shocked layers where the gas 
has cooled down sufficently to become molecular again. 
Alternatively, it could be due to ambient material put into 
motion by the passage of the atomic jet or to molecular 
gas in the outer layers of a disk wind. 



3. Derived physical parameters along the HHl jet 

Ratios of different lines can be used to constrain physical 
parameters, and their variation, along the jet, as described 
in some detail in Appendix A. However, to apply a diag- 
nostic procedure that combines lines that are relatively 
far apart in wavelength, as in this case, we need to deter- 
mine accurately the extinction towards each knot, and to 
correct the various line fluxes accordingly. To this aim, we 
have used the [Fell] lines at 1.32 and 1.64/im, that arise 
from the same upper level. The reasons for our choice, in- 
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Line id. 



A" 



F 



A F 



Line id. 



A° 



[Oi] ^D2- 
[Oi] ^D2- 
Ha 

[Nil] 

5/2 " 
3/2 



[Sii] 
[Sill 



3/2 
3/2 



[Fell] a^g/2-a^Fg/2 
[Fell] a\]r/2-a^Fj/2 
[Can] ^D,^2 
[Can] 2i?3/2 
[Nill] 2^7/2 -^i?5/2 
[Nill] 2^5/2- ^ 

[Fell] a^9/2 



-^'S'l/2 
'^'S'l/2 



^3/2 



a-^F 



[Fell] a4p5/2 



a'^7/2 
a'^3/2 



7/2 



[Fell] a4Pi/2 
[Fell] a^P3/2- ^'-^5/2 

[Cm] 6i59/2 -6^5/2 
[Cm] 6^7/2 -6^5/2 
[Cm] 6i55/2 -'^5/2 

[Cm] ^Dy2 -'^5/2 



3 /2 

[Fell] a4p5/2 
Can 2p°/2 - 



^5/2 



9/2 



[Fell] a% 



3/2 



3/2 

a"i^7/2 

[Fell] a''Pi/2~a%f2 
[Fell] a4p5/2 - a4p,/2 
[Fell] a^P3/2-a%^2 
[Fell] a4Pi/2 - 0^^3/2 

[Cl] li?2 -'Pi 
[Cl] li52 -=^^2 
[Sll]2 p3/2 -2 2^3/2 



6300 
6364 
6562 
6583 
6716 
6731 
7155 
7172 
7291 
7324 
7378 
7412 
7453 
7638 
7665 
7687 
8000 
8125 
8229 
8308 
8542 
8617 
8662 
8892 
9034 
9052 
9227 
9268 
9824 
9850 
1.029 



1.2410" 
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17 
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■15 


4.3810" 


17 


2.8610" 


■14 


5.8310" 
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1.8710" 
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17 
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17 
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16 
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17 
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17 
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17 
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15 
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16 
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16 
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16 
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15 
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16 


7.7510" 
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16 
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16 
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erg cm" 
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[Sll]2p3/2-2p5/2 
[Sll]2Pi/2-2p3/2 
[Nl]2P3/2,l/2 ~^D5/2 
[Nl]2P3/2,i/2 ~^D3/2 
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[Fell] a^D,^2 
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a^D^/2 
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1.189 
1.192 
1.248 
1.249 
1.258 



3/2 

- a'^5/2 

- a'^1/2 

- a^Dy2 

- a'^D^ 

4l 



3/2 



[Fell] a4p>3/2 
[Fell] a^Dy2 
[Fell] a^D3,2 
[Fell] a'^P'7/2 
[Fell] a^Dy2 
[Fell] a^Dy2-ci^F9/2 
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[Fell] a^Dy2-a^Fz/2 
[Fell] a^Dy2-a%/2 
H2 1-0 S(7) 
H2 1-0 S(6) 
H2 1-0 S(3) 
H2 1-0 S(2) 
H2 2-1 S(3) 
H2 1-0 S(l) 
H2 1-0 S(0) 
H2 2-1 S(l) 
H2 1-0 Q(l) 
H2 1-0 Q(2) 
H2 1-0 Q(3) 



1.280 
1.295 
1.299 
1.321 
1.329 
1.534 
1.601 
1.645 
1.665 
1.678 
1.748 
1.788 
1.958 
2.034 
2.074 
2.122 
2.224 
2.249 
2.407 
2.414 
2.424 
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Air wavelengths in A below 1/zm and vacuum wavelengths in micron above. 



stead of the more commonly used ratio [Feii]1.64/1.25^m, 
and the procedure followed, are described in Appendix B. 
The derived Ay values are listed in Table 2. There is a 
sharp decrease in the extinction, from Ay ~8 to 1.5 mag 
from the inner to the outer knots. Since for knots A to E 
the individual Ay values fluctuate more than the associ- 
ated errors, an average value was taken, with an uncer- 
tainty of 0.5 mag given by the dispersion of the different 
determinations. A visual extintion around 1.5 mag in these 
outer knots is in agreement with the E(B-V) value of 0.4 
found by Solf et al. (1988) towards the HHl bow shock, if 
a standard R=3 value is taken. 

As described in Appendix A, the selected optical tran- 
sitions from [S 11] , [O i] and [N 11] were used to determine 
the electron density, electron temperature, and ionization 



fraction in the knots adopting the BE99 technique, which 
in turn allow us to derive the total density nH=ne/xe ■ 

Given the limited spectral resolution of the EF0SC2 
configuration, the [S ll] doublet at 6716, 6731 A was not 
resolved sufficiently in all knots to allow us to separately 
measure the line fluxes. In these cases, we used data taken 
during a previous run with the same instrument but with 
higher spectral resolution, that allowed us to separate 
the doublet (Medves, Bacciotti & Eisloffel in prep.). We 
checked that the proper motion of the knots during the 
time interval between the two runs was negligible. 

The results of our analysis of the optical lines are 
shown in Tab. 3 and Fig.^ where the parameters are plot- 
ted as a function of the distance from VLAl. Moving out- 
wards, we measure a decrease in the electron density (from 
~10^ cm~^ in knots L-I to ~1.5102 cm"^ in knots C-B), 
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temperature (from 21000 K to 12000 K), as well as in the 



HH1 jet 



total density (from 1.210* cm"'* to -210'^ cm~ 



Such 



behaviour is often observed in HH jets (BE99). The ioni- 
sation fraction x^, on the contrary, is observed to slightly 
increase along the jet. This may be due to the fact that 
material ejected earlier had a higher degree of ionization 
and lower electron density. In fact, for the conditions found 
in stellar jets the ionization fraction is not in equilibrium 
with the gas thermal conditions, since the time scale for 
hydrogen recombination, which is proportional to l/rie, 
is generally longer than the travel time along the bright 
section of the jet. Therefore, the ionization fraction is es- 
sentially 'frozen' into the gas, and may trace past events 
at large distances. Alternatively, the gas may become re- 
ionized along the flow if, for example, the jet collides with 
an external clump (BE99) and an increase in ionization is 
induced by the UV photons created in the shock (Molinari 
& Noriega-Crespo 2002). This event, however, seems un- 
likely in our case, because the resultant shock should pro- 
duce an increase in the line fluxes which is not observed. 

The validity of our analysis can be checked by com- 
paring the fluxes of other observed lines against the values 
predicted assuming the derived physical parameters. Fig. 
[Slshows such a comparison for the N and S optical and NIR 
lines. The ratio [Sii](1.03Aim)/[Sii](6716,673lA) indicates 
the presence of higher temperature components. The same 
information is given by [NI](1.04/.tm)/[NII](6583,6548A) 
which, in addition, provides a check for the derived ioniza- 
tion fraction. In the BE99 diagnostics, in fact, the ioniza- 
tion fraction of hydrogen is found assuming that the O*' 
and N"*" populations are regulated by recombination and 
charge-exchange with H. Thus the observed [N I]/[N II] ra- 
tio can be used to verify this assumption. We see that for 
both [NI]/[NII] and [S ii](IR)/[S ii](opt) the agreement 
between the predicted and observed values is quite good 
for knots H and G while the predicted value is higher 
than the observed one for knots L-I. Since the same effect 
is seen in both line ratios, the problem does not arise from 
an overestimate of the hydrogen ionization fraction. Also, 
the presence of a temperature component higher than that 
derived by our analysis would produce a discrepancy be- 
tween the observed and predicted ratios even larger than 
the one we find. A better agreement between observed and 
predicted values in knots L-I can be obtained by assum- 
ing an Ay value smaller than estimated, i.e. ~6.5 mag 
instead of 8.3 mag. The agreement between observed and 
predicted ratio would also improve if the temperature in 
knots L-I is lower by about 60% with respect to the esti- 
mated value. 

The electron density and temperature can be measured 
in a completely independent way from the ratios of spe- 
cific [Fe ii] lines (see Fig. . Our results are also reported 
in Tab. 3 and plotted in Fig. 0] In Fig. |3 we see that the 
[Fe ii] model fits well all the observed [Fe ii] transitions in 
knot G. In fact the majority of the line ratios are repro- 
duced within about 30-40%. Likewise, the n,, determined 
from optical and NIR [Fe ii] lines decrease with distance 
from VLAl, however values derived using [Fell] are al- 
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Fig. 4. Derived physical parameters along the HHl jet. 
Open circles refer to parameters estimated using diagnos- 
tics based on [Fe ii] lines, while filled triangles refer to val- 
ues determined through the BE99 technique. Error sources 
are described in Table 2. 
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indicate the line ratios observed in the HHl jet knots and 
dereddened assuming the Ay values in Table 2. 
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Fig. 7. Dereddened [Fe ii] intensity ratios with respect to 
the 1.64/im line in knot G (open circles) against a model 
with re=11500 K and 71^=510^ cm^^ (filled triangles). 
Intensities are normalized to the 1.64/im line. 



from the [S ii] lines. Such a result has already been found 
in the analysis of other HH objects (Nisini et al. 2002), 
and indicates that the Fe emission originates from regions 
either denser and/or more ionized than those giving rise to 
the optical lines. The electron temperatures derived from 
the [Fell] line analysis are consistent with those derived 
from the optical diagnostics if we consider the associated 
uncertainty. The average values of Te([Fe ii]), however, are 
systematically lower than the Te(opt) values. In Fig.|31we 
see that such a behaviour is expected because of the dif- 
ferent spatial distribution of Fe and optical lines in the 
cooling region of shocks. In fact, Fe lines originate in a 
wide region of the post-shocked gas having, on average, a 
lower temperature and higher density than the zone pro- 
ducing the optical lines. 

3.1. Additional constraints on the excitation and 
physical structure 

As highlighted in Appendix A, several other line ratios can 
further constrain jet physical parameters. Our assumption 
of purely coUisionally excited gas can be checked using 
line ratios sensitive to fluorescence pumping effects. The 
dereddened [Niii]A7412/A7378 ratio is one of them, with 
expected values, for rie in the range 10^ to 10^ cm~^, 
between 0.04 and ~0.11 for pure collisions, and ~0.11 
and 0.4 if fluorescence excitation is included (Bautista et 
al. 1996). We obtain for this ratio a value ranging from 
0.04(±0.02) in knot H to 0.14(±0.06) in knot L-I, indicat- 



ing that the excitation is mostly coUisional. An additional 
probe of the absence of fluorescence pumping is given by 
the [Caii]A7291/A7324 ratio, which is expected to be 1.5 
for coUisional excitation (Hartigan et al. 2004). Our dered- 
dened ratios are always consistent with this value within 
the errors. This implies that the permitted Can transi- 
tions at 8542 and 8662A, detected in the inner knots E 
to L-I, are also excited by collisions. In turn, the ratio 
CaiiA8540/[Caii]A7291 can be used as diagnostic of high 
electron density as shown in Fig. |H1 This diagram has 
been constructed adopting a five-levels statistical equilib- 
rium code that uses radiative transition rates from NIST, 
and coUisional rates from Mendoza (1983) and Chidichimo 
(1981). This ratio is almost independent of the tempera- 
ture, while it starts to increase above ~0.1 only for elec- 
tron densities larger than ^10^ cm^^. The observed values 
thus imply rig ~10^ cm^^. 

High electron densities along the jet are also indicated 
by the [Feii]A7155/A8617 ratio, which, as described in 
Appendix A, is a tracer of regions at densities larger than 
~ 10^ cm~^. Figure shows the observed values against 
the predictions of Bautista & Pradhan (1998), for rie in 
the range 10^-6 10^ cm~^. Electron densities as high as 10^ 
cm~'^ seem to be reached in knot L-I, but a component 
at rie ~10^ cm^"^ persists further out along the jet. Note 
that these values are about one order of magnitude higher 
than those inferred from the NIR [Fell] lines. Values of 
electron density up to 10^ cm^^ have also been inferred 
in the inner and denser regions of few T Tauri stars by 
Hartigan et al. (2004). 
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Table 2. HHl jet: physical parameters 











Diagnostics from 0/S/N lines 


Diagnostics from Fe lines 


Knot 




Av ± AAt 

V V 


He ± Ane'^ 




hh ± Ann'^ 


Te ± ATe^ 


ne ± Ane'^ 


Te ± ATe'^ 




// 


mag 


10^ cm-3 




10^ cm-3 


K 


10^ cm-3 


K 


L-I 


1.2-3.8 


8.3±0.5 


2.8±0.1 


0.02±0.01 


101.1±40.6 


20 500±1900 


11±1 




H 


3.8-5.5 


2.9±0.5 


2.6±0.1 


0.06±0.01 


45.4±7.9 


13000±1100 


8±1 


9200t?r 


G 


5.5-8.5 


2.0±0.5 


2.1±0.1 


0.05±0.G1 


38.1±6.1 


11200±900 


5± 1 


10 5001?;^°° 


F 


8.5-11 


2.0±0.5 


1.2±0.1 


0.05±0.01 


26.4±6.2 


11700±1000 


4±1 


98001?^^° 


E 


11-12.4 


1.3±0.5 


0.9±0.1 


0.05±0.01 


17.7±5.1 


12800±1200 


30 ± 20 


8200^1^0 


D 


12.4 13.8 


1.3 ±0.5 


0.7±0.1 


0.08±0.G2 


8.2±3.2 


13 200±1400 




8300t}°S^ 


C 


13.8 15.6 


1.3±0.5 


0.6±0.1 


0.15±0.03 


4.1±1.9 


15100±2000 






B 


15.6-19.1 


1.3± 0.5 


0.7±0.1 


0.13±0.03 


5.4±2.1 


14 200±1800 




85001?^™ 


A 


19.1-26.2 


1.3±0.5 


0.6±0.1 


0.28±0.05 


2.0±0.8 


20300±2300 







"Distance from VLAl in arcsec 

''Visual extinction measured from the ratio [FeII]1.64/1.32/im (see Appendix 2 for details). 

■^uh is calculated as Xe/ug. Errors are estimated using both flux uncertainties and error propagation from uncertainties 
in the Ay determination and in solar abundance values reported in literature (Asplund et al. 2004). 
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Fig. 8. Calculated CaiiA8542/[Caii]A7291 ratio, as a 
function of the electron density and temperature. The 
solid lines indicate the range of values observed in the 
HHl jet. 

Finally, we point out the non-detection of the auroral 
optical [O ll] transitions in any of the knots. The 7331 A 
line pair are never detected, while the doublet at 732lA 
could be blended with [Caii]7324A. Since, however, the 
ratio [Caii]A7291/A7324 is always close to the value ex- 
pected from theory, the contribution of [O ii]A7321 should 
be negligible. This, in turn, is independent confirmation 
that the hydrogen ionization fraction, which is strongly 
tied to the 0+70° fraction by charge exchange, is low 
everywhere along the jet. In knot G, for example, the up- 
per limit [Oii]A7331/[Oi]A6300 -810-^ derived at the 3a 
level, implies an upper limit for the hydrogen ionization 
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Fig. 9. Observed values of the ratio [Fe ii]A 7155/A 8617, 
which is a tracer of high density regions. The two horizon- 
tal lines indicate the value expected for densities of 10^ 
and 6 10^ cm~'^ in a coUisionally excited gas at r=10000 
K (Bautista & Pradhan 1998). 

fraction Xe^O.l (Hartigan et al. 2004), which is consistent 
with Xe=0.05 derived by us. 

4. Mass flux in the jet 

The knowledge of the physical conditions along the jet 
beam allows us to derive the mass flux carried by the 
HHl jet, which is a fundamental quantity regulating the 
dynamics of the accretion/ejection process in protostars. 

For the determination of M we exploit two different 
procedures. In the first method, we directly estimate the 
mass flux in each knot combining the total density along 
the jet beam with values of jet radius, rj, and velocities, vj 
available in the literature (M = tuh nn xnVjVj, where 
niH is the proton mass, and /x=1.24 is the average atomic 
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weight). This method is not affected by uncertainties in 
the distance of the object and reddening estimate. On the 
one hand it assumes that the knot is uniformely filled at 
the density derived from the diagnostics, giving an upper 
limit to M. On the other hand, such an effect is partially 
compensated for by the presence of regions at even higher 
densities in the beams than those traced by the [S ii] lines. 

The mass flux derived with this method is given in 
Table 3 and plotted in Fig^] For rj, we have taken half 
of the FWHM of the [S ii] intensity profile in HST images 
measured by Reipurth et al. (2000), which increases from 
0.1 to 1 arcsec along the jet. Jet velocities are derived from 
proper motion studies by Bally et al. (2002), assuming an 
inclination angle of 10° (see Table 3). The derived mass 
flux ranges between 1 and 3 10~^ Mq yr~^ and is fairly 
constant along the jet within the uncertanties. 

The second method to derive the mass flux uses the 
observed luminosities of selected lines such as [S ii] , [O I] 
and [Fell]. These lines are optically thin, thus their lumi- 
nosity is proportional to the mass of the emitting gas (see 
e.g. Hartigan et al. 1984). We have: 

M = fimn X {uhV) X vt/k (1) 
uhV^ Lilme) (^/i z/ A, /, ^ (2) 

where Ai, fi are the radiative rate and fractional popula- 
tion of the upper level of the considered transition, X'^/X 
is the ionization fraction of the species X having total 
abundance with respect to hydrogen vt and It 

are the velocity and length of the knot, projected perpen- 
dicular to the line of sight. 

This method is affected by uncertainties in absolute 
calibrations, extinction, and distance, but does implic- 
itly take into account the volume filling factor. In prac- 
tice, the latter is given by the ratio of M values derived 
from the two different methods. We have computed M 
from the luminosity of [Sii]AA 6716,6731, [OI]A6300, and 
[Fe ii]1.64/im lines, assuming all Fe is in gaseous form 
(Tab. 4, Fig lTUI) . Since, as it will be discussed in section 
4.4, up to ~70% of Fe atoms can still be locked into grains, 
the derived M([Feii]) is actually a lower limit to the real 
value. 

Fig. 1101 illustrates the M values obtained with the dif- 
ferent methods. We note that the values derived from line 
luminosities are lower than those calculated from the jet 
density and radius, and fall off along the jet, indicating 
a decrease of filling factor. In fact, measuring the volume 
filling factor (//) by comparing the observed luminosity 
with that expected from the considered volume, one finds 
that // varies from ~1 to ~0.06 from knot L-I to knots 
A-B (Table 3). This behaviour is not unexpected, as the 
furthest knots have a larger volume filled with more dif- 
fuse gas, while the knots closer to the star are more dense 
and compact. 

In Fig. ^1 we note that the mass flux derived from 
[Fell] line luminosity is larger than M([Sii],[Oi]), at 



least from knot H outwards. This suggests that the 
higher values, given by the [Fe ii] diagnostic, come 
from gas with higher average density than in the [S ii] 
line emitting region, rather than from gas with a higher 
ionization fraction. Thus the [Fe ll] luminosity is probably 
a better tracer of the mass flux than the [O i] or [S ii] line 
luminosity, since the former traces a larger portion of the 
total mass flowing through a single knot. 

M derived from line luminosities, which takes into ac- 
count the filling factor, decreases with distance from the 
source. A constant mass flux along the jet is expected if 
the flow is stationary, since the decrease in the total den- 
sity should be compensated for by the larger jet radius. 
A mass flux decrease can be due to a change in the jet 
physical conditions: if the gas rccombines and the tem- 
perature decreases, the bulk of the flowing mass may be 
mainly in neutral atomic or molecular form and thus not 
be traceable through its ionic emission. To gauge the mass 
flux associated with the molecular components we have 
determined M(H2) using the 2.12/xm line luminosity. The 
temperature of the warm H2 gas is 2000-2500 K, as de- 
rived from the ratios of the dereddened line fluxes. Then, 
the total Nh2 column density is measured by comparing 
the intrinsic (per unit mass) and observed flux of the 1- 
S(l) 2.12/im line. The II2 tangential velocity has been 
assumed equal to that derived from optical lines, follow- 
ing Noriega-Crespo et al. (1997), who measured II2 knot 
proper motions comparable to those obtained from atomic 
emission lines. 

Our values of Af(Il2 ) range between 210^^" and 
2 10~^ Mqut^^, and are therefore about 2-3 orders of 
magnitude lower than those derived from the atomic lines. 
Moreover, M(H2) does not increase with distance. Thus it 
is unlikely that the decrease of the mass flux derived for 
the atomic/ionic component is due to atomic gas turned 
into molecular material. Other contributions, however, 
such as atomic gas at low temperature, not emitting in 
the optical and near IR, may become important to the 
bulk of mass flow. In this regard, we note that strong 
[OI]63/im line emission, which traces shocked atomic gas 
at temperatures below 5000 K, has been detected by 
ISO in a 70" beam centred on VLAl (Giannini et al. 
2001, Mohnari & Noriega-Crespo 2001). Similar 63^m 
emission has also been found towards HHl and HH2, 
suggesting that such emission is not excited in the source 
circumstellar envelope, but derives from both the shocked 
gas and the diffuse FDR of the hosting cloud. This 
latter contribution has been estimated to account for 
about 50% of the total [01] luminosity. Excluding the 
FDR contribution and estimating the mass flux from 
the [OI]63//m luminosity (HoUenbach 1985), we derive 
M{[Ol]63nm) - 310^6 Mgyr-^ Such a large value 
suggests that a significant part of the mass flux may be 
due to a warm atomic component. FIR observations at 
resolution higher than that available with ISO, may be 
obtained in the near future with SOFIA and Herschel, 
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Table 3. Mass flux in the jet and fiUing factors 



knot 




Tjet 


l\/f ['riTT r - 

ivi [nn, r-j. 


-1 


Ml\^ ttI W 


f ffK tt1"1« 
//IP IIJJ 


MtlVn ttI W 


-ivi© yr 




km s ^ 


ouL csec 


J-Vi© yr 


J-vj-o yr 




J-Vi© yr 


T T 

Ij-i 


oUU 


U.i 


i.o iU 


-7 


0/11 n^7 




1 n^7 




H 


300 


0.15 


1.6 10" 


-7 


3.6 10"* 


0.2 


7.1 10-8 


4.710-1° 


G 


305 


0.2 


2.410" 


-7 


6.910-^ 


0.3 


1.010-^ 


1.910-9 


F 


300 


0.3 


3.110" 


-7 


3.710-s 


0.1 


5.210-8 


8.910-1" 


E 


281 


0.3 


2.310" 


-7 


1.810-8 


0.08 


4.4 10-8 


1.310-9 


D 


300 


0.35 


1.510" 


-7 


1.510-8 


0.1 


3.810-8 


1.810-9 


C 


296 


0.4 


9.710" 


-8 


1.110-8 


0.1 


2.910-8 


1.610-9 


B 


317 


0.5 


1.710" 


-7 


7.510-9 


0.04 


2.410-8 


2.110-1° 


A 


311 


0.5 


7.610" 


-8 


5.110-9 


0.07 


2.810-9 


1.110-9 



" taken from the proper motion study by Bally et al. (2002). Total velocities {vjet) have been derived from these 
values assuming an inclination angle of 10°. 

^ taken from the FWHM intensity transverse profiles in HST images of Reipurth et al. (2000) 

Af (riH, I'j.Vj) = ^niH nH x Trr^gj Vjet, where uh is the total density reported in Table 2. Volume filling factors equal 
to unity are assumed. 

M measured from the luminosities of the [Sii]AA 6716,6731 lines and [Feii]1.64/im line. 
^ Volume filling factor estimated from the ratio between the observed and predicted line luminosities (see text). They 
are also equal to the ratio between M and M{nYi,rj,Vj) derived from fine luminosities. 

M measured from the luminosity of the H2 2.12/im line, assuming the H2 temperature estimated from the ratio of 
the different H2 observed lines. 



and will allow us to gain a more complete picture. 

Finally, we note that M(H2) in the HHl jet is on aver- 
age lower by one or two orders of magnitude with respect 
to the values derived by Davis et al. (2001, 2003) for other 
small-scale H2 jets from embedded sources. These authors 
also find that M(H2) and M(Fen) are roughly comparable 
for the majority of the objects in their sample. In contrast, 
our estimates for HHl show its this jet, despite its origin 
from an embedded source, presents characteristics similar 
to jets from more evolved YSOs. In such sources the molec- 
ular component contributes little to the jet excitation and 
dynamics in comparison to the atomic component. 

We also note that the molecular outflow associated 
with the VLAl source is weak and less energetic than 
other outflows from Class sources (Moro-Martm et al. 
1999). The momentum flux deduced by our observations 
of the atomic jet (^3 10-''' M© km yr~^) is comparable 
to that estimated for the associated CO outflow (Correia 
et al. 1997). All thi evidence suggests that around VLAl 
the density of molecular material is low and thus that the 
dynamics of the flow is mainly determined by the atomic 
gas. 

5. Abundances of refractory elements 

The observed wavelength range includes several lines from 
refractory species. In addition to Fe, we detected transi- 
tions from C, Ca, Ni and Cr. These species are expected 



to be largely depleted in neutral or molecular interstellar 
regions, while their gas phase abundance is enhanced in 
shocked regions where the dust grains have been at least 
partially destroyed either by sputtering or by photoevap- 
oration (e.g. Draine 2003). Therefore, the determination 
of gas phase abundances can provide constraints on dust 
properties and structure, on dust destruction mechanisms 
and on how such mechanisms depend on physical parame- 
ters such as jet velocity and excitation conditions. In this 
regard we note that the different species are expected to 
follow selective patterns for their erosion from the dust 
grains, and thus the fraction of their gas-phase abundance 
can be different and give further indications of dust grain 
reprocessing. 

To derive quantitatively the gas-phase elemental abun- 
dance, however, one has to rely on comparing emission of 
refractory and non-refractory species that originate in the 
same region. In this way, the difference between predicted 
ratios of lines, calculated with the determined values of 
Te, Ue and Xe, and the observed ratios gives indications 
on the depletion of the refractory species with respect 
to the assumed solar abundance. For Fe this procedure 
is quite complicated, and contrasting results have been 
reported for HH objects (see e.g. Bohm & Matt 2001, 
Beck-Winchatz et al. 1996, Nisini et al. 2002). The fact 
that the Fe spectrum is consistent with physical condi- 
tions different from those inferred from the optical lines 
(e.g. [Sii]A6716 and [Oi]A6300) gives strong limitations to 
a diagnostics based on a comparison between such lines. 
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This is demonstrated in Figure 1111 which is a 
plot of the hnc ratios [Fe ll]A8617/[S ii]AA6716,6731 and 
[Feii]A8617/[Oi]A6300, assuming a Fe gas-phase abun- 
dance equal to solar (310~^, Grevesse & Sauval 1998). 
The [Feii]A8617 line has been chosen because its spatial 
distribution is similar to the distribution of the [S ii] and 
[O i] lines, as shown in Figure 3. However, FigurelTTlshows 
that the predicted ratios, assuming for all lines the phys- 
ical conditions derived from the BE99 analysis, underes- 
timate the observed values and, in fact, would suggest Fe 
gas phase abundances larger than solar. 

We have then considered the ratio 
[Feii]1.25/um/[Pii]1.18/im (Figure [HI). Phosphorus 
is a non-refractory species and its 1.18/im line has about 
the same excitation conditions as the [Feii]1.25/im line. As 
shown in Oliva et al. (2001), the ratio [Feii]1.25/[Pii]1.18 
should be about [Fe/P]/2. Assuming solar abundances 
([Fe/P]~ 112) we find [Fe]ga,/[Fe],oiar -0.7 for knots G 
and F and ~0.2-0.3 for the internal knots L-I and H, in 
line with values derived in Nisini et al. (2002) for other 
HH objects. Finally, a Fe gas-phase abundance lower than 
solar is also inferred from a comparison of [Fe ii] lines and 
Pa/3 as in Nisini et al. (2002). In knot G, taking a Pa/3 
upper limit of 5 lO'^^ ergs-^cm-^, Te=11000 K and 
j:e=0.05, we derive a lower limit of [Fe]gas=0.1 [FeJ^oiar 
from the plot in Fig. 13 of Nisini et al. (2002). 

The gas phase abundance of C and Ca is more eas- 
ily derived from a comparison with the [S ii] optical lines 
assuming the same excitation conditions (see Fig. 3). 
Figure 1131 shows the comparison between the observed 
and predicted ratios [Ci]AA9824,9850/[Sii]AA6716,6731 
and [Caii]AA7291,7324/[Sii]AA6716,6731. 

While Ca can be assumed to be completely ionized, 
the fraction C'^/Ctot has been computed with an ion- 
ization balance between collisional and charge-exchange 
ionization, and direct/dielectronic recombination (rates 
from Standi et al. (1998) and Landini & Monsignori Fossi 
(1990)). Figure US shows that the observed [Ci]/[Sii] 
ratio decreases towards the source, in contrast to what 
can be expected theoretically. In knots L-I, however, the 
higher temperature causes C to become almost completely 
collisionally ionized and the theoretical ratio is also re- 
duced. The discrepancy between observed and expected 
[Ci]/[Sii] ratios can be interpreted as due to a decrease 
in the C gas-phase abundance. This behaviour can be cor- 
related to the increase of the total density from about 5 10^ 
to 1.2 10^ cm^^ going from knots A-C to knots L-I, since 
the efficiency of depletion increases with the density. On 
the other hand, the higher temperature of knot L-I may 
cause dust evaporation of C to compete against depletion. 

In Fig. El we also show that the observed [Caii]/[Sii] 
ratio is smaller than predicted for knots F to L-I, al- 
though Ca depletion is lower than for C. To check if 
selective depletion is present among other species, we 
have plotted in Figure US the [Niii]A7378/[Feii]A8617 
and [Cr ii]A8000/[Feii]A8617 ratios observed in the various 
knots. There is evidence that the [Niii]A7378/[Feii]A8617 



ratio decreases with distance. Since Ni+ and Fe"*" coex- 
ist, and the two lines have the same excitation conditions 
(Bautista et al. 1996), we interpret the observed decrease 
as a change in the [Ni]/[Fe] gas abundance ratio. In this 
respect it is interesting to note that Beck-Winchatz et al. 
(1996) found a Fe abundance close to solar in the HHl 
bow shock, but here the Ni/S ratio was 10 times the solar 
value. 

The presented analysis, therefore, shows that a certain 
degree of depletion of the refractory elements is present, 
an effect which is apparently selective, being higher for 
Fe and C. This shows that dust grains are not completely 
destroyed in stellar jets, and that the number of refrac- 
tory atoms released in the gas phase is an inverse function 
of the total density. In fact, models for dust reprocessing 
in shocks show that shattering processes, and grain-grain 
collisions, are able to destroy only a few percent of large 
grains for shock velocities Vg ^ 10 kms~^ (Draine 2003). 
At the same time, sputtering in fast shocks produces com- 
plete dust destruction for Vg > 100km s~^ (Jones 2000). 
Since shocks in stellar jets, like the HHl jet, have veloc- 
ities no larger than 40-50kms^^, we would expect that a 
large fraction of dust grains survive shocks and contribute 
to the observed partial depletion of refractory elements. 
Moreover, our results also show that multiple low velocity 
shock events, such as those occurring in protostellar jets, 
are not 100% efficient in completing the dust destruction 
process. 
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Fig. 11. Ratios of [Feii]8617A line with respect to 
[Sii]AA 6716,6731 and [Oi]A6300 in the various knots. 
Crosses: measured de-reddened values. Filled circles: pre- 
dicted values assuming and derived from the BE99 
analysis. Fe and S have been assumed completely ionized 
once, while the O+ZO" fraction has been calculated using 
the ionization model for this species. 
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Fig. 13. Dereddencd ratios of lines from refractory 
species (crosses): (a) [Ci]AA9824,9850/[Sii]AA 6716,6730, 
(b) [Caii]AA7291,7324/[Sii]AA 6716,6730, (c) 

[Nin]A7378/[Fcii]A8617, (d) [Crn]A8000/[Fcii]A8617. In 
panels (a) and (b), filled circles are the predicted values 
assuming Ca, C and S have solar abundances, C all 
neutral and Ca, S all singly ionized. Open circles in panel 
(a) refer to predictions where the C^/C"*" abundance 
ratio has been computed through ionization equilibrium 
analysis. 



6. Summary 

We have presented a new combined optical/NIR analysis 
of flux calibrated spectra of the HHl jet. The simultane- 



ous use of several diagnostic lines has proven to be very 
powerful in probing the density and temperature stratifi- 
cation inside the jet knots, and has allowed us to obtain a 
very detailed picture of the jet physical structure and exci- 
tation even from moderate spectral and spatial resolution 
observations. 

Using different line ratios, we have derived visual ex- 
tinction, electron density, temperature and ionization frac- 
tion as a function of distance from the driving source. The 
electron density from the [Sii] doublet decreases from 
310^ cm-3 in the L-I knots (at d<4" from VLAl) to 
~600 cm"'' in the A-D knots far from the star. Conversely, 
Ue values between 10^ and 4 10'^ cm"'' are derived in the 
same knots from [Fell] line IR ratios. Finally, the use of 
very high density probes, such as [Fe ii] A7155/8617 and 
CaiiA8542/[Caii]A7291, supports the presence of jet com- 
ponents at even higher electron densities, ranging from 10^ 
to 10^ cm^''. Such a stratification presumably originates 
in the different layers of a post-shocked gas. This is reason- 
able since the post-shock cooling distance is expected to 
range between ~10"'^^ and 10^^ cm (Hartigan et al. 1987, 
see also Fig. 1), i.e. smaller than our spatial resolution. 
Our conclusion is also supported by the temperature val- 
ues inferred from the [Fell] optical/NIR ratios, which are 
lower in comparison to those derived from using optical 
line diagnostics. Fe emission therefore traces post-shocked 
regions located at larger distances from the shock front 
than [S ii] lines, where the compression is higher and the 
temperature is declining. 

The ionization fraction, derived using the BE99 diag- 
nostics, is rather low (xe ~0. 02-0.1) which is also sup- 
ported by the non-detection of the [Oil] auroral lines at 
7321 and 733lA. Several line ratios have been also used to 
exclude the presence of fluorescence excitation even in the 
knots closest to the source. This implicitly validates the 
diagnostic procedures we use, which assumes that the gas 
radiative properties are regulated by coUisional excitation. 

If we assiime that the ionization fraction is almost con- 
stant throughout most of the post-shocked line emitting 
region, as the profile in Fig. 1 would suggest, we deduce, in 
the innermost jet region, a total density as high as 5 10^ 
cm"''. Moreover, the range of total densities probed by 
our analysis implies that in each knot a compression fac- 
tor C = ripostshock/fiproshock of at Icast 10^ is reached. For 
comparison, values of C lower than '--^100 are derived from 
dissociative shocks with velocity ranging from 30 to 100 
kms~^ (Hartigan et al. 1994). 

The mass fiux in the jet has been determined for each 
knot using two methods. In the first M is calculated com- 
bining the estimated Uh values with the jet radius and 
velocity taken from the literature {M{nH,rj,Vj)). In the 
second we used a comparison of observed and predicted 
total luminosity of the [Sii], [Oi], and [Fell] lines. In the 
jet section closest to the source (i.e. knots L-I) both meth- 
ods suggest a value of ~210~'' M0yr~^. Conversely, as 
we move away from the driving source we find discrep- 
ancies between the results from the two methods. The 
M{nH, fj,Vj) value represents an upper limit to the mass 
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flux, since one makes the assumption that the jet is com- 
pletely filled by gas at the estimated density. Conversely, 
M([SII],[OI]) calculated from line luminosities, takes into 
account the proper filling factor, but underestimates the 
real M value, since in this case one does not include the 
higher density gas regions. M([FeII]) is intermediate be- 
tween these two values, and is probably a better tracer of 
the global mass flux through the jet. 

Af ([Fell] is found to decrease with the distance from 
the driving source, from 2 10~^ in knots L-I down to ~3- 
5 10~^ M0 yr^^ in knot A. Such a decrease may indicate 
that in the more distant regions of the jet an atomic com- 
ponent at T < 5000 K (as traced by the [OI]63/zm emission 
observed by ISO) becomes important and contributes to 
the bulk of the mass flow. The mass flux due to the warm 
molecular H2 gas has also been estimated and it turns out 
to be always smaller than the other values by one or two 
orders of magnitude. 

Finally, the presented observations suggest that the re- 
fractory elements are still partially depleted with respect 
to the assumed solar abundance by fractions which varies 
from less than 10% in the low density regions up to '^90% 
for C, ~70% for Fe and ~50% for Ca in the inner and 
densest knots. This is a result expected in relatively low- 
velocity shocks, since the complete destruction of dust by 
shattering processes is predicted only for shock velocities 
larger than 50-100 kms~^. It also indicates that the ac- 
tion of several low- velocity shock events during the jet life- 
time is not effective in producing a complete destruction 
of dust grains. While we found that there is evidence of 
depletion depending on gas density, an analysis performed 
on a larger number of jets is needed to properly determine 
the dependence of the gas-phase abundances of refractory 
elements on jet excitation conditions. In addition, the ob- 
servation of a selective pattern for the depletion of differ- 
ent elements needs to be investigated further since it may 
provide important information on dust grain composition 
and properties. 

Appendix A: Diagnostic capabilities of the 
observed lines 

A.l. Oxygen, Nitrogen and Sulphur lines 

The forbidden lines from O, N, S in the optical range 
are powerful diagnostic tools for stellar jets, as dis- 
cussed in BE99, Dougados et al. (2000), Hartigan et 
al. (2004). In particular, the ratio of [Sii]A6716 and 
A6731 is sensitive to electron densities in the range 50- 
210^ cm^-^, this latter value being the critical den- 
sity for coUisional excitation. Furthermore, BE99 have 
developed a method to derive the hydrogen ioniza- 
tion fraction Xe and the electron temperature Te from 
the observed ratios [O i]AA6300,6363/[Nii]AA6583,6548 
and [Oi]AA6300,6363/[Sii]AA6716,6731. The technique is 
based on the fact that in typical conditions for the beams 
of stellar jets the hydrogen ionization fraction is primarily 
determined by charge-exchange with oxygen and nitro- 



gen atoms, while collisional ionization of H and photoion- 
ization do not make a large contribution until the shock 
velocity exceeds 100 kms^^. Such a velocity, however, is 
only reached in the large bow-shocks at the head of the 
jets, while those in the beam are generally not larger than 
30-40 kms~^. An additional assumption is that S is all 
singly ionized. Such an assumption is motivated by the 
large LP. of S+ (50 eV) and by the non-detection of the 
[Siii] and [Si] lines at 9533A and 1.082/xm respectively, 
in any of our spectra. More details can be found in BE99 
and Bacciotti (2002). The jets analysed with this tech- 
nique in previous studies typically have 0.01< Xe <0.6 
and 7000< Te <210^. When using the BE technique a 
set of elemental abundances has to be assumed. In this 
study we have adopted solar abundances from Grevesse & 
Sauval (1988). 

[S 11] and [N i] lines are also observed in the NIR part of 
the spectrum. They consist of four closely spaced lines for 
each species ([Sii] at 1.029, 1.032, 1.034, and 1.037/Ltm 
and [Ni] at 1.0400, 1.0401, 1.04100, and 1.04104/im) 
having excitation energies higher than the optical lines 
(^35 000 and 40000 K, respectively). Thus, in principle, 
the [Sii]1.03//m lines can be combined with the optical 
[Sii] hues to give an independent estimate of Tg. On the 
other hand, the [Ni] lines at 1.04^m can be used with the 
[N 11] optical lines to determine the ionization fraction. 

A. 2. Iron lines 

Most of the observed [Fell] lines originate from transi- 
tions among the first 16 fine structure levels. As illus- 
trated in Nisini et al. (2002), electron densities in the 
range lO'^-lO^ cm~^ can be derived from line ratios of IR 
lines, 1.64/1.60/im and 1.64/1.53/im in particular, while 
the ratios 1.64/1.25/im and 1.64/1.32/i,m can be used 
to determine the visual extinction (see Appendix B). In 
the optical part of the spectrum, we also observe tran- 
sitions coming from levels at higher excitation tempera- 
tures. Consequently, the ratios between selected optical 
and IR [Fe 11] lines are a good diagnostic tool for the mea- 
surement of the electronic temperature in the range 5 000- 
20 000 K (Nisini et al. 2002, Pesenti et al. 2003). 

The combination of the IR and optical spectrum of 
ionized Fe can therefore be used effectively to give an in- 
dependent diagnostic of temperature and density, which, 
relying on ratios from lines of a single species, is not af- 
fected by any assumption of abundances. For our analysis, 
we use a statistical equilibrium model which considers the 
first 16 [Fell] levels, the collisional strengths from Zhang 
& Pradhan (1995) and the transition probabilities from 
Nussbaumer & Storey (1988). An example of a diagnostic 
diagram for such an analysis is shown in Fig. Similar 
diagrams are also presented in Nisini et al. (2002) and 
Pesenti et al. (2003). 

We have also detected, in some of the knots, the lines 
at 7155, 7172, and 7453A which connect the doublet 
0^^9/2,7/2 with the a'^F levels. As shown in Bautista & 
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Pradhan (1998) and Hartigan et al. (2004), the ratio be- 
tween the 7155A and 8617A transitions is sensitive to 
electron densities in the range ^^10^-10^ cni~^, sig- 
nalling the presence of high density conditions. 

A. 3. Other atomic/ionic lines 

Other prominent atomic/ionic lines observed in the spec- 
tra include transitions from the refractory species [C I] , 
[Call], [CrII] and [Nil!]. 

The [C I]AA9824,9850 lines can be excited both by col- 
lisions with neutrals and electrons in a partially ionized 
medium, or by recombination (Escalante & Victor 1990). 
In this latter case, other lines from the recombination cas- 
cade, such as the 1.07 and the 1.17/im transitions, which 
are expected to be only a factor 1-10 weaker than the 
9840A lines, should also be detected (Walmsley et al. 
2000). Their non-detection in any of our spectra would 
imply that the [C I] lines are collisionally excited in the 
HHl jet. The presence of strong [CI] lines also implies 
low excitation conditions, and in fact the ratio [C I] /Ha is 
a very good indicator of excitation in HH objects. 

We detect [Call] lines originating from transitions 
among the first five Ca+ levels. In particular, we detect 
the two forbidden lines at 7291 and 7324A connecting the 
^^5/2.3/2 levels with the ground state, as well as the per- 
mitted lines at 8542 and 8662A, connecting the ^^3/2, 1/2 
levels with the ^£'5/2, 3/2- These latter two transitions can 
be collisionally excited, being only at 25 000 cm~^ from 
the ground state or, alternatively, photoexcited through 
the Call K and H lines. Thus, the ratios among the ob- 
served lines can be used as a diagnostic to see if fluorescent 
excitation is important. 

Ni+ is expected to coexist with Fe"*", given their sim- 
ilar ionization potential (7.9 and 7.6 eV, respectively). It 
has been found to be overabundant with respect to the 
solar value both in HH objects and in other gaseous neb- 
ulae such as supernova remnants (Beck-Winchatz et al. 
1996, Bautista et al. 1996). In the latter, the observed 
strong [Niil] emission has been explained by the contri- 
bution of fluorescent excitation in a low density environ- 
ment. Indeed, the ratio between the optical lines, at 7412 
and 7378A can be used as a discriminator between colli- 
sional and fluorescent excitation (Bautista et al. 1996). 

Finally, in the spectra of the brightest knots we also 
detected the two [Pii] transitions at 1.148 and 1.189/xm. 
Such lines have excitation properties very similar to the 
[Fell] lines but P is a non-refractory element. This makes 
the IR ratio [Feii]/[Pii] a good tool to estimate Fe deple- 
tion (Oliva et al. 2001). 

Appendix B: Determination of Ay from [Feii] 
lines 

Visual extinction can be estimated from pairs of lines orig- 
inating from the same upper level, providing the emission 
is optically thin. In this case, their theoretical intensity ra- 
tio depends only on frequency and transition probabilities, 



and not on physical conditions in their emission region. 
There are several [Fell] lines which have this property, 
and the ratio [Feii]1.64/1.25^m (hereafter 1.64/1.25) is 
most commonly used for this purpose, as these lines are 
the brightest of the near IR [Fell] spectrum (Gredel et 
al. 1994, Nisini et al. 2002). Among detected [Fell] lines, 
the bright transition at 1.32^m also comes from the same 
upper level as the 1.64 and 1.25/im lines, and thus the 
ratio [Feii]1.64/1.32/^m (hereafter 1.64/1.32) gives an in- 
dependent check of reddening. In Fig. IB. II we plot the ob- 
served 1.64/1.25 and 1.64/1.32 ratios in the different jet 
positions. This plot shows the reddening variations along 
the jet irrespective of the absolute value of the visual ex- 
tinction, indicating that the reddening is almost constant 
from knots A to H and sharply increases in knots L-I. This 
is consistent with Reipurth et al. (2000) who estimated 
that the visual extinction in the inner portion of the jet 
is larger, by ~ 4 mag, than that in the external knots. To 
convert the observed ratios to Ay values, we adopt the 
radiative transition probabilities given by Nussbaumer & 
Storey (1988), which imply intrinsic ratios 1.64/1.25=0.73 
and 1.64/1.32=2.81. We further use an interpolation of 
the near IR Rieke & Lebofsky (1985) extinction law of 
the form Ax/ Ay = 0A2 \{pm)-^ ''^ (Draine 1988). In 
Figure IB. II we also plot the Ay values derived from the 
two line ratios. Note that the extinction estimated from 
the 1.64/1.25 ratio is always about a factor of two larger 
than the values estimated from the 1.64/1.32 ratio. 

Such differences can be critical when estimating the 
intrinsic ratio of lines far apart in wavelength. Thus, we 
need to investigate the possible causes of such a discrep- 
ancy in more detail. 

The observed discrepancy cannot originate from any 
systematic observational effect. The three lines (1.64, 1.25, 
and 1.32/im) are all observed with the same grism and thus 
are not affected by any problem of intercalibration. They 
are detected at large S/N, implying that any contamina- 
tion due to OH sky line residuals (which could affect the 
1.64/xm line) is negligible. We then checked for observa- 
tions of these three lines in other objects. In Nisini et al. 
(2002) the Ay values estimated from the 1.64/1.25 ratio 
in a number of HH objects were always higher than the 
value derived from optical observations. Moreover we have 
verified that extinction values derived from the 1.64/1.32 
ratio are more consistent with the optical determinations. 
The same effect, i.e. Ay (1.64/1.25) > Ay (1.64/1.32), can 
be found in various gaseous nebulae where these lines 
have been detected with sufRcent accuracy (e.g. Orion B, 
Walmsley et al. 2000, and supernova remnants, Oliva et 
al. 1990). Finally, in the HHl bow shock Noriega- Crespo 
& Garnavich (1994) estimated an Av = 6.7 using the 
1.64/1.25 ratio, pointing out the discrepancy between this 
value and that measured through optical line ratios. 

The most likely explanation for such a discrepancy is 
probably uncertainty in the transition probabilities of the 
lines involved. More recent calculations of such probabili- 
ties, by Quinet et al. (1997), show discrepancies with re- 
spect to the Nussbaumer & Storey (1988) values as large 
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Fig. B.l. Observed variations along the HHl jet of the 
reddening sensitive [Fell] line ratios 1.64/1.25/xm (a) and 
1.64/1.32/xm (b). Panel (c) shows the derived Ay values 
from the two ratios assuming an interpolation of the Rieke 
& Lebofsky (1985) extinction law and transition probabil- 
ities from Nussbaumer & Storey (1988). Circles and trian- 
gles are the Ay values measured from the 1.64/1.25 and 
1.64/1.32 line ratios, respectively. 



as 20%. With the Quinet et al.(197) A values, the the- 
oretically expected ratios become 1.64/1.25=0.965 and 
1.64/1.32=3.67 ^. If applied to our data, such ratios pro- 
duce on average Ay smaller than the values derived with 
the NS88 data, but the discrepancy between the 1.64/1.25 
and 1.64/1.32 estimated values remains. In knot G, for 
example, we derive, with the Quinet et al. (1997) values, 
Ay(1.64/1.25)= 1.6 mag, but Av(1.64/1.32) < 0. In gen- 
eral, the yl\/(1.64/1.32) value estimated with the param- 
eters by Quinet et al. (1997) tends to be always negative, 
suggesting that the theoretical value is probably overesti- 
mated. 

Finally, we point out that the higher Ay values derived 
from the 1.64/1.25(NS88) ratio, turn out to be inconsis- 
tent with the ratios of other observed transitions irrespec- 
tive of the assumed physical conditions. For example the 
[Feiijl. 64/0. 862 ratio in knot G, if corrected for an Ay 



In a note added in proof, Quinet et al. (1995) evaluated the 
limitations of their Aij calculations done considering only three 
[Fe ll] electronic configurations. While in general the inclusion 
of the additional 3d^4d^ configuration does not change the ra- 
diative rates by more than few percent, for the explicit case 
of the a*F-a*D transitions (which includes the 1.64^m transi- 
tion) the correction is about 40%. If real, this would cause the 
1.64/1.25 theoretical ratio to become again very similar to the 
1.64/1.25 NS88 expected ratio. 



larger than 4, is outside the range of expected theoreti- 
cal ratios for any temperature (see Fig. 4). The same ap- 
plies to the ratio [S ii]0.67/1.03 (see also the discussion in 
Section 4.2). In view of the above discussion, we suggest 
that the 1.64/1.32 ratio from NS88 gives, for the moment, 
the most reliable determination of the Ay value. 
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